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Abstract 


An  asymptotic  analysis  is  given  of  the  compressible,  laminar 
boundary -layer  flow  of  a  dilute  gas-particle  mixture  over  a  semi-infinite 
flat  plate.  The  analysis  extends  existing  work  by  considering  more 
realistic  drag  and  heat-transfer  relations  than  those  provided  by  Stokes.  A 
more  general  vi scosi ty-temperature  expression  is  also  incorporated  into  the 
analysis.  The  solution  involves  a  series  expansion  in  terms  of  the  slip 
parameter  of  the  particles.  The  numerical  results,  including  the  zeroth  and 
first-order  approximations  for  the  gas  and  particle  phases,  are  presented 
for  the  two  limiting  regimes:  the  large-slip  limit  near  the  leading  edge 
and  the  small-slip  limit  far  downstream.  Significant  effects  on  the  flow 
produced  by  the  particles  with  Stokes'  and  non-Stokes'  relations  are  studied 
and  clarified.  The  effects  of  some  nondimensional  similarity  parameters, 
such  as  the  Reynolds,  Prandtl  and  Eckert  numbers,  on  the  two-phase 
boundary-1 ayer  flow  are  discussed. 
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coefficient  in  Eqs.  (5.4)  and  (5.20) 

coefficient  in  Eq.  (5.4) 

coefficient  in  Eqs.  (5.5)  and  (5.21) 

coefficient  in  Eq.  (5.5) 
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total  heat  transfer  per  unit  volume  to  gas  from  particles 
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1.0  INTRODUCTION 


1 . 1  Motivation  for  the  Present  Study 

Gas  and  solid-particle  flows  are  encountered  in  many  different  fields. 
Typical  examples  occurring  in  nature  are  dust  storms,  forest-fire  smoke  and 
the  dispersion  of  solid  pollutants  in  the  atmosphere.  Many  processes  in 
industry  utilize  gas-particle  flows,  such  as  transportation  of  pulverized 
materials  in  pneumatic  conveyors,  separation  and  classification  of  particles 
in  cyclone  or  other  separators,  fluidization  in  chemical  reactors,  and 
combustion  of  powdered  fuels  in  combustion  chambers.  In  addition,  gas  flows 
with  suspended  solid  particles  have  various  applications  in  science  and 
engineering,  for  example,  satellite  drag,  ablation,  MHD  generators,  solid 
propellant  rockets,  1  aser-Doppl er  anemometry  and  blast  waves  moving  over  the 
Earth's  surface. 

For  some  applications  in  pipe  or  nozzle  flows  and  Hows  over  bodies, 
the  behaviour  of  such  two-phase  flows  at  a  solid  surface  is  extremely 
important.  Hence,  it  is  necessary  to  study  boundary-layer  flows  of  a 
gas-particle  mixture.  From  solutions  of  gas  boundary-layer  equations,  it  is 
possible  to  determine  the  effects  of  solid  particles  on  the  boundary-layer 
characteristics,  say,  shear  stress,  heat  transfer  and  boundary-layer 
growth . 

The  problem  considered  in  this  report  is  the  laminar  boundary-layer 
flow  over  a  semi -i nf i ni te  flat  plate  in  a  compressible  gas  containing 
uniform,  spherical  solid  particles.  This  study  provides  basic  physical 
insight  into  the  flow  of  such  a  two-phase  system,  even  though  the  solution 
is  for  a  relatively  straightforward  problem.  Moreover,  as  a  parallel  study, 
the  asymptotic  solution  can  be  used  to  compare  with  finite-difference 
solutions  and  to  verify  independently  the  correctness  of  the 
finite-difference  scheme  [1]. 


1.2  Previous  Work 


Several  authors  have  worked  on  the  problems  of  two-phase  boundary-1 ayer 
flows.  Most  of  these  analyses  were  based  on  the  assumption  of  an 
incompressible  fluid  [2-16],  Singleton  [17]  first  treated  the  case  of  a 
compressible  dusty-gas  boundary-layer  flow.  He  derived  the  governing 
equations  and  obtained  asymptotic  solutions  for  two  limiting  regimes:  the 
large-slip  regime  near  the  leading  edge  and  the  small-slip  regime  far 
downstream.  However,  he  assumed  Stokes'  relation  for  the  drag  force  and 
heat  transfer,  which  is  valid  only  for  the  case  where  the  particle  slip 
Reynolds  number  is  of  order  unity.  He  developed  his  governing  equations 
assuming  that  the  gas  viscosity-temperature  relation  has  the  special  form  of 

u*/u*  =  /T*/T*,  and  gave  his  solutions  for  the  case  where  the  Prandtl  and 
Eckert  numbers  of  the  gas  are  equal  to  unity. 


^  k  .1 


AmJI 


•/vVv/;A 
,  A  A  A  A  A 


A*  A 


l 


1.3  Present  Study 

The  present  analysis  will  extend  Singleton's  analysis  to  the  more 
general  problem  of  compressible  laminar  dusty-gas  boundary-layer  flows  over 
a  semi-infinite  flat  plate.  It  will  present  the  basic  equations  under 
conditions  that  the  drag  and  heat  transfer  between  the  two  phases  may  have 
different  relevant  forms  instead  of  Stokes'  relation  and  that  the  power 
index  in  the  expression  for  the  viscosity  coefficient  can  have  arbitrary 
values  from  0.5  to  1.0.  The  paper  will  give  the  numerical  results  in  the 
two  limiting  regions  at  several  values  of  the  Prandtl  number,  Eckert  number, 
Reynolds  number  and  the  viscosity  power  in.ie<. 


1.4  Basic  Assumptions 

The  basic  assumptions  are  as  follows: 

(1)  The  gas  is  perfect.  The  specific  heats  of  the  gas  are  constant.  The 
Prandtl  number  of  the  gas  is  constant.  The  viscosity  and  heat 
conductivity  of  the  gas  have  a  power-law  relation  with  the  gas 
temperature. 

(2)  The  solid  particles  are  rigid  spheres  of  uniform  size.  The  number 
density  of  particles  is  sufficiently  high  to  treat  the  particle  phase 
as  a  continuum.  However,  the  particles  are  also  sufficiently  dilute  to 
consider  them  as  non-interacting. 

(3)  The  particles  have  no  random  motions  and  therefore  the  particle  phase 
does  not  contribute  to  the  static  pressure  of  the  two-phase  system. 

(4)  The  volume  fraction  of  the  particle  phase  is  assumed  as  negligible. 
This  implies  that  the  coefficient  of  viscosity  for  the  gas-particle 
mixture  can  be  taken  as  the  viscosity  of  the  gas  phase  alone. 

(5)  The  specific  heat  of  the  particle  material  is  constant.  Its  thermal 
conductivity  is  much  larger  than  that  of  the  gas  and  hence  the 
temperature  inside  each  particle  can  be  assumed  uniform. 

(6)  There  is  no  radiative  heat  transfer  from  one  particle  to  another. 
There  is  no  chemical  reaction,  no  coagulation,  no  phase  change  in  the 
two-phase  system.  There  is  no  particle  deposition  on  the  surface  of 
the  plate. 

(7)  Only  the  processes  of  drag  and  heat  transfer  couple  the  particles  to 
the  gas.  The  drag  coefficient  and  the  Nusselt  number  for  a  single 
sphere  in  a  viscous  flow  are  assumed  valid  for  the  particle  cloud. 
Other  force  interaction  terms,  such  as  lift,  buoyancy  and  gravity,  are 
neglected. 
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(8)  There  is  no  dry  friction  as  the  particles  slide  along  the  wall.  The 
slowing  down  of  the  particle  motion  is  only  due  to  the  gas  whose 
velocity  decreases  to  zero  at  the  wall. 

(9)  The  two-phase  flow  is  steady.  The  flow  Reynolds  number  is  sufficiently 
high  so  that  a  laminar  boundary-layer  forms  on  the  surface  of  the  flat 
plate,  but  lower  than  a  critical  value  so  that  no  transition  to 
turbulence  occurs. 

(10)  The  usual  boundary-layer  assumptions  are  still  valid  for  the  two-phase 
system  and  consequently  the  variation  of  pressure  across  the  boundary 
layer  can  be  neglected.  In  addition,  for  the  flat-plate  problem,  there 
is  no  pressure  gradient  in  the  external  flow.  The  particle  phase  and 
gas  phase  in  the  external  flow  are  in  equilibrium. 

2.0  GOVERNING  EQUATIONS  AND  BOUNDARY  CONDITIONS 
2.1  Governing  Equations 

Let  x*  and  y*  be  the  distance  along  and  normal  to  the  wall, 
respectively.  The  origin  is  fixed  at  the  leading  edge  of  the  plate.  The 
geometry  of  the  problem  is  sketched  in  Fig.  1.  The  conservation  equations 
for  steady  two-dimensional  laminar  boundary -layer  flows  of  a  compressible 
gas-particle  mixture  over  a  semi-infinite  flat  plate  are  as  follows: 

For  the  gas  phase: 

Continuity: 


® _  p*u*  +  - —  p*v*  =  0 

dx*  dy* 


Momentum: 


p*  (u*  +  v*  ^u*  )  =  — _  (p*  ^u*  )  +  Dx 

dx*  by*  by*  by* 


Energy: 


State: 


p*  =  p*R*T* 


(2.1) 

(2.2) 

(2.3) 

(2.4) 


I  For  the  particle  phase: 

i 

i 

j  Continuity: 


!  _  p*  ui  +  - c£  v*  =  0 

5x*  HP  P  ay*  T  P 

(2.5) 

x-momentum: 

.  au*  au* 

p*(u*  — E.  +  v*  — E.)  =  -Dx 

Pv  P  3x*  P  by*1 

(2.6) 

i 

y-momentum: 

p*(u*  +  V*  =  -Dy 

Pp  p  ax*  p  ay*  y 

(2.7) 

|  Energy: 

j 

(2.8) 

In  Eqs.  (2.1)— (2.8) ,  the  interaction  terms 
particles  can  be  expressed  by  (see  Appendix): 

between  the  gas  and 

! 

u*  -  u* 

Dx  -  $  D 

V 

(2.9) 

y*  _  y  ^ 

d*  *  <$  -V- 0 

(2.10) 

<1>  =  (u*  -  u*)Dx  +  (v*  -  v*)D y 


(2.11) 


where 


TS  -  T* 


Q  =  p*  c*  P  Nu 


(2.12) 
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Tf  = 
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(2.16) 


Here,  D  represents  the  real  drag  coefficient  Cq  normalized  by  the 
Stokesian  drag  coefficient  Cq0  and  Nu  is  the  Nusselt  number  based  on  the 
particle  diameter.  According  to  the  assumption  (7),  D  and  Nu  determine  the 
gas-particle  interaction.  When  the  interaction  law  between  the  gas  and 
particles  just  has  the  Stokes  form,  D  =  1.0  and  Nu  =  2.0.  It  is  well  known 
that  the  Stokes  relation  is  valid  only  for  small  slip  Reynolds  number  of 
order  unity.  In  general,  D  and  Nu  are  both  functions  of  the  Reynolds  number 
and  Prandtl  number.  From  the  definition  of  slip  Reynolds  number  and  Prandtl 
number. 


Re, 


p*/(ug  -  u*)2  +  (vg  -  v*)2  d* 


(2.17) 
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(2.18) 


The  local  equilibrium-time  parameters,  x*  and  xf,  are  a  measure  of  the 
relaxation  process.  For  example,  the  velocity  equilibrium  time  x*  is  the 
time  elapsed  for  a  particle  to  reduce  its  relative  velocity  to  e- 1  of  its 
original  value  if  the  force  accelerating  (or  decelerating)  the  particle 
toward  the  gas  velocity  is  given  by  the  Stokes  drag.  These  two  local 
parameters  are  functions  of  the  local  gas  temeprature  since  the  viscosity 
coefficient  p*  and  the  heat  conductivity  k*  are  functions  of  the  gas 
temperature.  It  is  convenient  to  introduce  an  equilibrium  length  \*,  which 
is  based  on  the  freestream  parameters: 


.  Ps 


H*  2 


18  p* 


(2.19) 


The  two-phase  relaxation  process  takes  place  throughout  the  equilibrium 
length.  Therefore,  it  is  reasonable  to  choose  X*  as  the  characteristic 
length  of  the  dusty-gas  boundary-layer  problem.  Then  the  equilibrium-time 
parameters  can  be  expressed  in  the  form 


(2.21) 


•>  C*  X*  |i* 

tf  =  —  Pr  —  *  - —  *  — 

T  2  c*  u*  p* 

Thus,  the  basic  boundary-layer  equations  (2.1)-(2.8)  become 

— _  p*u*  +  — _  p*  v*  =  0 

dx*  ay* 


P*  (u*  OP*  +  V*  — *  )  =  (p*  +  fV*(u„  -  u*)  —  -H- 

p  1  3v*  ftv*  1  P  X*  .1* 


ay*  ay 


p*(u*  HI  +  v*  HI)  =  i_  ±-  (n*  H 1)  +  j£  (&£) 
ax*  ay*  Pr  ay*  3y*  c*  ay* 


v-p  ^ao  r^oo 

+  JL  p*(T*  -  T*)  °J*»1  Nu 
3Pr  PV  P  X*  M* 


p*  =  p*R*T* 


J! p*  u*  +  Pi*  v*  =  0 

3^  PP  P  3y*  T>  P 


au*  auS  ,  ^  ,  u$, 

U*  — B.  +  V*  — P-  =  -  u*  -  u*)  —  —  D 

p  ax*  P  ay*  p  x*  p* 


U*  +  V*  =  -(v*  -  V*)  —  H.  D 

P  ax*  P  ay*  p  x*  u* 


U*  3.  ♦  V*  =  -  _5L  (T*  -  T*)  "1  M*  Nu 
p  ax*  P  ay*  3Pr  P  X*  p* 


(2.22) 


(2.23) 


(2.24) 


(2.25) 


(2.26) 


(2.27) 


(2.28) 


(2.29) 


my  JsmmmmmwTmn 


where  a  is  the  ratio  of  specific  heats  of  the  two  particles: 


a  =  _£ 

r* 


(2.30) 


For  boundary-layer  flows,  the  normal  component  of  velocity  is  usually  a 
small  quantity.  It  means  that  the  contribution  of  the  normal  velocity  is 
often  neglected  compared  with  the  tangential  velocity.  Then,  the  expression 
for  the  slip  Reynolds  number,  Eq.  (2.17),  becomes 


Res  = 


p*]u*  -  u*|d* 


(2.31) 


and  the  gas  energy  equation  (2.24)  reduces  to 


p*(u*  511  *  v*  511)  »  ]_  (p»  511)  +  JC  (55l)2 

H  av*  ;  Pr  v  r*  ftv* 


ay*  Pr  5y 


by*  c*  ay^ 


+  (u*  -  u*)2  ~  —  D  +  (T*  -  T*)  —  Ji*  Nu 

c*  P  X*  p*  3Pr  P  X*  p* 


(2.32) 


In  order  to  obtain  a  closed  set  of  equations,  it  is  required  to  specify 
the  expression  for  p*(T).  From  the  standpoint  of  kinetic  theory  of  gases, 
the  most  exact  relation  for  the  viscosity  of  a  perfect  gas  is  Sutherland's 
form  [18]: 


=  15-121  £)3/2 

p*  T*  +  S*  'o 


(2.33) 


where  S*  is  the  Sutherland  constant.  This  form  of  Sutherland's  relation, 
however,  is  not  suitable  for  the  series-expansion  method  which  is  used  in 
the  present  analysis.  The  other  viscosity  form  used  in  many  analytical 
solutions  to  the  boundary-layer  equations  is  of  power  form,  which  is  written 


4  -  (£f 


(2.34) 


where  w  is  the  power  index  which  lies  between  0.5  and  1.0.  This  relation 
for  the  gas  viscosity  is  readily  applied  to  the  series-expansion  method,  as 
shown  later. 


m 


m 


m 


From  the  basic  assumption  (10),  the  pressure  is  constant  throughout  the 
boundary  layer: 


p*  =  constant 


(2.35) 


With  this  condition  (2.35)  and  the  gas  state  equation  (2.25),  the  gas 
density  can  be  expressed  in  terms  of  the  gas  temperature: 


*  T* 

P*  _  CD 

£  T7 


(2.36) 


Substituting  Eq.  (2.36)  into  the  basic  equations  (2.22),  (2.23), 

(2.26)-(2.29)  and  (2.32)  with  the  power  relation  for  the  gas  viscosity 
(2.34),  the  following  equations  are  obtained: 


W*7!  +  *1)  ,  u*  HI  +  v* 

ax*  by*  dx* 


(2.37) 


u*  sm  ♦  V*  M  =  v*„(H)  »_  [(H)"  &£]  .  (U*  -  u*)  ^  (Hf*1  D 

•V  .  4.  **v  X  C"VTO.  '  -v  -±-  *-  Wj,  7  ^  i  J  *  p  '  4.  'T*  ' 


T*  ay*  T*,  ay 


T* ' 


(2.38) 


u*  .511  +  V*  .511  =  —  (H)  1_  r (11)“  -5111  +  1®  (II)^1  (M) 

ax*  ay*  Pr  T*  ay*  UT*  ay*  c*  T*  0y* 

*  4s  <“*p  -  u*>2 1  ^  I  <T5  -  T*>  £  Nu  (2-39> 


3Pr  p* 


X*  T* 

'vOD  1  00 


.5 p*  u?  +  — —  p?  v*  =  0 

fix*  P  P  fiv*  T  P 


u*M+v*^=-(u*-u*)^(Il)% 
P  ax*  P  ay*  P  X*  ^t*; 


(2.40) 


(2.41) 


S'  .*  v  '.-.v  v  v  .v  -• : 


w; 


•-  i  .A.  i* -A  T  - 


(2.42) 


Sv* 

U*  — ^  + 

P  ax* 


9v* 

v*  = 
P  ay* 


.(v*  - 
1  P 


V*) 


^  [lift 

X*  T* 


dT* 

u  *  — £  + 
P  dx* 


u*  ,T*.W 


V*  m  =  -  _JL  (T*  -  T*)  If  (Ill 
P  ay*  3Pr  P  X*  lT* j 


Nu 


(2.43) 


Under  the  conditions  of  w  -  0.5,  D  =  1.0  and  Nu  =  2.0,  the  equations 
(2.37)-(2.43)  reduce  to  those  derived  by  Singleton  [17]. 

Physically,  the  boundary-layer  flow-field  of  a  two-phase  mixture  can  be 
divided  into  three  distinct  regions  (see  Fig.  1).  These  regions  are  divided 
according  to  the  nondimensi onal  slip  parameter  x*/X*  as  follows:  the 
large-slip  region  (x*/X*  <<  1),  the  moderate-slip  region  (x*/X*  ~1),  and 
the  small-slip  region  (x*/X*  >>  1).  Following  Singleton  [17],  a 
smal 1 -parameter  expansion  method  was  used  to  solve  the  boundary-layer 
equations  for  dusty  gases.  Only  the  asymptotic  solutions  in  the  two 
limiting  regimes  can  be  obtained  by  this  perturbation  technique:  the 

large-slip  approximation  for  the  near  leading-edge  solution  and  the 
small-slip  approximation  for  the  far-downstream  solution,  respectively. 
Clearly,  the  large-slip  regime  is  characterized  by  a  frozen  flow  where  the 
gas  and  the  particles  move  independently,  while  the  small-slip  regime  is 
characterized  by  an  equilibrium  flow  where  the  gas  and  the  particles  move 
together  (see  Fig.  2) . 

2.2  Boundary  Conditions 

The  boundary  conditions  for  the  gas  phase  are: 

(1)  At  the  wall,  there  is  neither  slip  in  velocity  nor  jump  in 
temperature: 

u* ( x* ,  0)  =  0,  v* ( x* ,  0)  =  0,  T*( x*,  0)  =  T*  (2.44) 

(2)  As  y*  approaches  infinity,  the  flow  parameters  must  match  those  in  the 
external  flow  or  the  freestream: 


u*(x*,  ®)  =  u*,  T*(x*,  »)  =  T*  (2.45) 


The  boundary  conditions  for  the  particle  phase  are: 
(1)  At  the  wall,  there  is  no  mass  transfer 
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V*(x*.  0}  =  0 


(2.46) 


(2)  As  y*  approaches  infinity,  the  flow  parameters  must  match  their 
freestream  values: 


u*(x*,  •)  =  u*^. 


T£(x*,  -)  -  T*m, 


p£(x*.  -)  =  p*^  (2.47) 


Since  the  particles  and  gas  are  assumed  to  be  in  equilibrium  in  the  external 
flow,  the  freestream  parameters  for  the  particle  phase  can  be  readily 
determined  as 


T*  =  T* 

'  p  1  <»* 


=  3p* 


(2.48) 


where  p  is  the  mass  loading  ratio  of  the  particles.  Otherwise,  the 
two-phase  external  flow  must  be  solved  first  in  order  to  obtain  the  outer 
boundary  conditions  for  the  particle  phase  if  the  particles  are  not  in 
equilibrium  with  the  gas  in  the  freestream. 


3.0  LARGE-SUP  APPROXIMATION 


3.1  Transformation  of  Boundary-Layer  Equations 


a 

>5 


na 


u  - 


X 

x* 

=  -L_ ,  T)  = 

V5, 

/  u»  * 

/  2v*  x*  y 

(3.3) 

V  = 

/  2X*  V*. 

/  V*  u*, 

T  =  — ,  p  =  £ 

T*  P* 

(3.4) 

f  =  nu  - 

v’  fp  -  J 

/  1  **„ 

2v*u*x*  P 

(3.5) 

PP 

„  PL= 

PL  pp~’ 

T* 

TP  =  y6 

oo 

(3.6) 

- 

P  -  4- 
A 


(3.7) 


With  Eqs.  (3.1)-(3.3)  and  (3.5)-(3.6),  the  particle  velocities  can  be 
expressed  as 


<  =  u*[i-  -J>] 


Pp  aii 


__  [1_  &  ♦  f  -  IL  !!e)] 

v  Pp  3x  2x  2x  dn 


(3.8) 


(3.9) 


Substituting  the  above  expressions,  the  basic  equations  (2.37)-(2.39) 
and  (2.41-(2.43)  are  transformed  into  the  following  form: 


t r au_  +  u_-l_^L)  =  u  5L.I _ ar 

ax  2x  2x  an  ax  2x  an 


£_au_T__a_  rjw  Su )  +  (°  P 

2x  an  2x  an  an  an 


Ppu  )D 


(3.10) 


(3.11) 


+  y_  -  1_  at  =  —L  £_  (T03  *L)  +  T"^)*  +  pEc  &  -  o  U  Y 

dx  2x  2x  3n  2xPr  do  an  2x  an  o-  an  ™ 


*  Wr  PPT“(TP  '  T,NU 


Pp  8r, 


^Le  [«  -  jl  ^i)  -  ^e  (^e.  -  jl  1^2.)] 

dn  P  ax  drj  2x  qv2  an  ox  2x  an 


(^E  +  !e  -  : L  !!e)(p  ^Le  -  ^E  ^E.) 

ox  2x  2x  an  P  ar|2  an  an 


-  vtV  PpU)D 


.p  .liA.i  £lz  *  jl  !!ej 

p  an  ax2  x  ax  4x2  2x  ax  an  4x2  an 


f  dPp  ^Pp  _  Yp  ^Pp  U^P  +  ^JP  _  n  ^p  ) 

an  ax  ax  an  2x  an  ax  2x  2x  an 

+  p_(^£  +  !e](^!e  -  jl  £!i)  =  pD^rw(ile-  +  Ie  -  jl  ^!e  +  jl  P  u 

P  dx  2x  3xdn  2x  3^2  P  3x  2x  2x  dn  2x  P 


^!e  !Ie  -  (!!e  +  !e)  = .  _e_  P  tw(t  -  t)nu 

an  ax  an  ax  2x  3Pr  P  P 


where  Ec  is  the  gas  Eckert  number  based  on  the  freestream  temperature 


Ec  =  — =  (Y  -  U^2 

r  *  T* 

°° 


3.2  Expansion  in  Terms  of  the  Slip  Parameter  (x*A*) 


The  following  expansion  in  terms  of  the  slip  parameter  (x*/\*)  is 


with  the  boundary  conditions 


f(0)(0)  =  0;  u(0)(0)  =  0;  u(0)H  =  1;  T^°> (0)  =  Tw;  t(°)(«)  =  1  (3.24) 

fp0)  (0)  =  0;  fj0)‘(»)  -  1;  Tj0)(-)  -  1;  pj°)(»)  =  l  (3.25) 

In  fact,  the  equations  for  the  particle  phase,  (3.21 )-( 3.23) ,  are  algebraic 
expressions  which  represent  the  zeroth-order  solution  for  the  particles. 
The  boundary  conditions  (3.25)  were  already  used  in  Eqs.  (3.21 )-(3.23) . 

The  first-order  problem  is 


f(ir  _  rl«V  f(D  _  3u(0  t  (2u(0)  ,  Tmv_  f(0),  tUT  .  f<°>  T(i) '  -  n 

T(0)  T<0)  T«*>  T<°>  (3.26) 

U(1)"  *  (_r<£L  +  -  2  "(0)  ..(1)  -  [(<^1)  f(0)u(0),i 

T(0)wfi  t(°)  _tn\ufr2 


T(0} 


T(0) 


. ,,  t<°)'u(°)'  ]T(i)  t  ..  U(Q)-  T»).  .  U(Q>-  ,m 
2  j(0)  f n\w*-l 


T(0) 


T(0) 


=  -20(1  -  u<°>  )D 


(3.27) 


t(D"  +  (Pr  _I 


(0) 


T(0) 


u>+- 1 


T(0) 


T(0) 


w+-2 


j(0) 


+  2Pr  ij^°- — ]T  ^ 1  ^  +  2EcPr  (u ^  'u ^  '  )  +  Pr  f^) 

_(0)wfl  /mi**! 


•(0) 


=  -20EcPr(l-u(°)  )2D  -  l  0(1  -  T<°)  )Nu 


(3.28) 


3f 


-  3^)'  +  3ft1)  =  -2T(u)U(t1u(0)  -  f(°))D 


(3.30) 


nTj1)'  -  T ( 0 ) T ( 0 )  -  l)Nu 

P  P  3Pr 


(3.31) 


where  the  normalized  drag  coefficient  D  and  the  Nusselt  number  Nu  are  given 
by  their  zeroth-order  approximations .  The  boundary  conditions  are: 

f(1)(0)  =  0;  uW(0)  =  0;  u(1)M  =  0;  T^^O)  =  0;  T^H  =  0  (3.32) 

f^1}(0)  =  0;  f^1)'(«)  =  0;  Tp1)(“)  =  0;  p^M-O  (3.33) 

From  the  first-order  approximate  equations,  Eqs.  (3. 26)-(3. 31 ) ,  it  is 
seen  that  the  loading  ratio  of  the  particles  p  appears  oly  in  the  equations 
for  the  gas  phase,  i.e.,  Eqs.  (3.26)-(3. 28) .  In  addition,  it  is  possible 
to  get  a  more  general  solution  which  is  suitable  for  any  value  of  the 
loading  ratio,  by  introducing  the  new  variables: 


F(1)(n)  =  f(1)(n)/(3 


(3.34) 


U(1)(r,)  =  u<1>(r1)/p 


(3.3b) 


(1)U)  =  T(1)(h)/P 


(3.36) 


Substituting  the  expressions  (3. 34)- (3. 36 )  into  Eqs.  (3. 26 )- (3. 28 ) ,  the 
first-order  equations  for  the  gas  phase  become 

u).  (1,  ,i)  ,o)  tJ(ur  ,0))_m  t  fW  U). 

T(0>  T(0)  J  T(0)  T(0) 


uto-  *  1 ) ■  -  2  -i(-°2-uU) 

po)"*1  i(0>  T(o)“tl 


Conti nued 


©(!)“  ♦  (pr  f<0)T(0>'  , 

T(0)u,fl  T^°)  T  (O)^2 

+  2Pr  Jj.(_L.]©(1) 

T(O)^1 


-2EcPr(U(0)  V1^  )  -  Pr  -I - ^  -  2EcPr(l-u(°)  )2D  -  -  (l-T^0)  )Nu 

(O)^1  3 

'  (3.39) 


Similarly,  the  boundary  conditions  (3.32)  take  the  form 


F(1)(0)  =  0;  l)( 1  ^  ( 0)  =  0;  U^H  =  0;  ©(^(O)  =  0;  ©(^H  =  ^  ^ 


3 . 3  First-Order  Problem  and  Interaction  Terms 

From  Eqs.  ( 3 . 18) -( 3.23)  and  (3.26)-(3.31) ,  it  is  found  that  the 
interaction  terms  between  the  gas  and  particles  appear  only  in  the 
first-order  approximation.  In  other  words,  the  zeroth-order  equations  can 
be  solved  without  knowing  the  interaction  relation  between  the  two  phases. 
However,  in  order  to  obtain  the  first-order  solution,  appropriate 

expressions  for  the  drag  and  heat  transfer  between  the  two  phases  should  be 
gi ven. 

As  mentioned  before,  when  solid  particles  move  through  a  gas  at  very 
low  relative  velocities,  that  is,  when  Re$  <  1,  the  Stokesian  form  can  be 

applied.  For  the  case  where  Stokes'  relation  applies,  then  D  =  1.0  and 

Nu  =  2.0.  When  the  relative,  or  slip,  velocities  between  the  particles  and 

the  gas  increase  to  a  higher  value,  Stokes'  relation  is  not  valid. 
Therefore,  a  relevant  form  for  the  interaction  should  be  assumed  for  the 
non-Stokesian  case.  In  the  case  of  larger  slip  velocities,  it  may  be 


reasonable 
[19,  20]: 


to  apply  the  followiny  dray  and 


heat-transfer  relations 


c0  =  0.48  +  28  Re"0*85 

(3.41) 

Nu  =  2.0  +  0.6  Pr 1/3  Res1/2 

(3.42) 

The  drag  coefficient  Cq  and  Nusselt  number  Nu  for  the  non-Stokesian 
case  in  this  analysis,  given  by  Eqs.  (3.41)  and  (3.42),  are  functions  of  the 
slip  Reynolds  number  Res  as  well  as  the  Prandtl  number  Pr.  Comparing  the 
non-Stokes  and  the  Stokes  relations  indicates  that  these  two  cases  are  in 
agreement  only  for  the  slip  Reynolds  number  of  order  unity  or  less.  As  the 
slip  Reynolds  number  becomes  larger  than  unity,  the  Stokes  relation 
underestimates  the  drag  and  heat-transfer  between  the  gas  and  the  particles. 
In  general,  the  non-Stokes  relation  agrees  with  the  standard  dray  curve  [21] 
much  better  than  the  Stokes  relation.  Correspondingly,  the  normalized  dray 
coefficient  is 


D  = 


+  1  Re,.0*15 
6  s 


(3.43) 


In  the  series-expansion  method,  the  slip  Reynolds  number  Re$  should  be 
expanded  just  as  the  other  quantities.  Neglecting  the  first-order  small 
quantities,  the  slip  Reynolds  number  can  be  expressed  as 


Re$  =  Rep  l-  -u- 


(0) 


j(0) 


OJ+1 


(3.44) 


where  Rep  is  the  particle  Reynolds  number  based  on  the  freestream  velocity 


Re 


P 


d* 


p 


★ 

ao 


(3.45) 


The  zeroth-order  approximations  for  D  and  Nu  can  be  obtained  by 
substituting  Eq.  (3.44)  into  Eqs.  (3.42)  and  (3.43).  Then  the  first-order 
equations  can  be  solved  numerically.  The  equations  for  the  yas  phase,  Eqs. 
(3.26)-(3.28) ,  consist  of  second-order,  ordinary-differential,  simultaneous 
equations  with  two-point  boundary  values,  similar  to  the  case  of  the 
zeroth-order  equations  for  the  gas  phase.  The  solution  to  the  equations  for 
the  particle  phase,  Eqs.  (3. 29)-(3. 31 ) ,  can  be  obtained  in  the  integral 
form: 
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v. 


vvjSP,‘ 


& 


V-V.-.o.s 


(3.46) 


Pp1)  *  -2 rj2  /  —  —  (x  -  f(°))Ddx  -  n 2  /  li-i—  (xu^0)  -  f(°))Ddx 

U 

a>  X  ®  X 

♦  /"  liil!  (xu<°>  -  f(0))Ddx 

00  X  2 


f]l)'  .  -3nJ  (xu<°)  -  f<°>)Ddx  ♦  (xu<°)  -  f«»)Ddx 

®  x4  ®  x2  (3.47) 


f(l)  .  V  /"  l!£l!  (xu(°>  -  f(°)  *>dx  V  n  /"  li—  (xu<°>  -  f<°>  »dx 

“X  “>  x2  (3.48) 


T^)  =  -  j  T^_  (T(0)  .  i  )Nu  dx  (3.49) 

3P  r  qo  ^  3 

The  quantity  f^1)  =  dfj^/drj,  given  by  (3.47),  can  be  used  to  give  the 

first-order  approximation  of  the  tangential  velocity  for  the  particle  phase, 
ui1'.  From  Eq.  (3.8),  the  x-component  of  particle  velocity  can  be  given  by 
the  derivative  of  the  transformation  function  fp(n): 

1  Sfn 

Up  =  —  — ^  (3.50) 

M  Pp  ^ 


In  addition,  in  the  large-slip  limit,  the  nondimensional  density  of  the 
particle  phase,  pp ,  is  of  order  unity  since  the  zeroth-order  solution  is 
py'  =  1,  i.e.,  Eq.  (3.23).  Substituting  Eqs.  (3.17),  (3.21)  and  (3.23), 
Eq.  (3.50)  yields  the  series-expansion  form  as 

up  =  1  +  x(f(D’  -  p^))  +  ...  (3.51) 


Clearly,  the  first  and  second  terms  in  Eq.  (3.51),  represent  respectively 
the  zeroth  and  first-order  approximation  of  the  tangential  particle  velocity 


4.0  SMALL-SLIP  APPROXIMATION 

4.1  Basic  Equations  in  Terms  of  Slip  Quantities 

For  the  small-slip  region,  it  is  convenient  to  employ  slip  quantities 
as  dependent  variables  since  they  are  small  quantities  of  first  order  with 
respect  to  the  slip  parameter  The  slip  quantities  are  defined  as 


u*  =  u*  -  u*. 


T*  =  T*  -  T* 
s  p 


(4.1) 


Putting  Eq.  (4.1)  into  the  basic  equations  (2.37)-(2 .43 )  and  making 
some  algebraic  manipulation  results  in 
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The  conditions  of  D  =  1.0  and  Nu  =  2.0  are  already  employed  in  Eqs. 
(4.2)-(4.8),  since  the  slip  velocities  are  always  small  quantities  in  the 
small-slip  region. 


4.2  Transformation  of  Basic  Equations 
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As  in  the  large-slip  limit,  substituting  Eqs.  (4.9)-(4.12)  into  Eqs 
(4.2)-(4.8),  the  dimensionless  basic  equations  can  be  obtained: 
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where  x  =  x*/X*  is  the  nondimensional  slip  parameter. 


Then,  the  small-slip  approximation  for  the  zeroth-order  problem  leads 
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with  the  boundary  conditions 


f<0)(0)  =  0;  u(0)(0)  =  0;  u(0)(»)  =  1;  T(0)(0)  =  Tw;  T  ^ 0  ^ ( “)  =  1 

(4.25) 


For  the  small-slip  approximation,  the  zeroth-order  solutions  of  the  particle 
velocity  and  temperature  are  the  same  as  those  for  the  gas  and  the 
first-order  solutions  are  given  by  the  first-order  slip  quantities.  The 
zeroth-order  density  for  the  particle  phase  is  given  by  Eq.  (4.24).  The 
boundary  condition  for  the  particle  phase  density 
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has  been  used  during  the  derivation  of  Eq.  (4.24). 

The  first-order  equations  for  f^,  u^^,  T^  and  p(^  are  given  by 


f(l)'  .  T(0) '  f(l)  =  .  T(0) 1  f(0)T(l)  +  2  T(l)  +  T(l)'  .  u(l) 

T(°)  t(0)2  T(0)  T(0)  (4 


u(l)"  +  f(Q)  +  u( 1) '  +  2u(0)  u(l) 

T(0)wfl  T(o) U)fl 

__  _p_  1  .  2  f (QMQU0)'  +  (hM)  f(^^  ^t(Q)  1 

Up  T(0)^l  2T(0^  7(0)“  T(0)Wfl 


+  f(°)3u(0),]  .  u(O)1  f(l)  .  _J_  f(0)u (0)' 
2T(0)2wfl  T(o)wfl  1+p  T(0)wfl 


(liil  +  T(°}q<1)) 
T(0)  * 


♦  .1^1(1)  ♦  (*1)  1M1t(1)  -  ^T(l)'  (4.28) 

T(0)2  T(0)wf2  T(0) 

T(1)h  ♦  (Pr  1  VP/?f(0)  +  Z^WV0)')  -L(Q)-1  +  [2Pr 

UP  ;  T ( o)  ^ 1  1+P 


Conti nued 


24 


V.v  V  .'.v  />. .-v 


U  'llflUR 


*T 


-  Pr  ^l  *  P/g)  +  1  t(Q)'  f(°)  -  T(0)<  2  ] 

1+P  T(0)  T(0)1_w  jCO)"4-1 


R: 

i 


s 

i> 


5 


=  pr  p/a  1  rnu(0^  T(0)l  .  f  (Q)u(Q)T(Q)  ’  +  3(2o>-l)Pr/  a-2  ft°)  2J 
1+P  /n\W  fn\U  4  /  n\ 


2_ln\<  2 


2i(°) 


2t(°) 


T(0) 


(0) 

W+-1 


9 Pr  f(0)u(Q)T(°)1  +  3Pr 2  1  +  p/  a  f(°)  T ^ 0 ^ '  +  3EcPr2  fWy0)'  ] 
4a  t(0)u>  4  a  1+  P  T(Q)2wfl  4a  T(0)w 


EcPr  JL  '“W* 

1+P 


i ‘ 2  _  Pr  i  +  p/q  tC°) '  f(i) 

2T(0)2^  1+P  y^)^1 

-  2EcPr>(°)  V1)'  ) 


Pr  P/g  f(Q)T(Q)'  t1) 

1+P  T(o)“  ^ 


(4.29) 


Lv 

LV 

M 

ty. 


fOlp'11'  *  (2u<°>  *  lEl  f  (»)  )p‘l) 

t(0) 


'Jl 

i-r 


sf 


■\v 


:S 


f (o)u(Q) 1  +  nu(Q)u(0)t  +  ^  f(o)u(o)i(Q) 1  +  U  f (Q) ^(Q) 1 

2  T(0 )Wf2  ytO)^3 


ylO)^1  T(0)Wfl 


+  i  [Pr  1  •*•  P/«  f  (°>T<0> '  +  EcPrfu^)’  )21  f(°)2  -  i£tl  nu(Q)  2T(0)‘ 

2  1+P  T(0)u>fl 


rfo)1^2  2  yfo)^2 


+  f '  T(l)  _  2  u^0)  T(l)  f(Q)  y(l) • 

T(0)3  T(o)2  t(°)2 


(4.30) 


with  the  boundary  conditions 
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The  first-order  problem  in  the  small-slip  limit  is  determined  by  the 
second-order,  ordi nary-di f ferenti al ,  simultaneous  equations  with  two-poinc 
boundary  values.  In  this  aspect,  it  is  similar  to  the  first-order  problem 
in  the  large-slip  limit.  But  the  solution  to  the  above  equations 
(4.27)-(4.30)  is  difficult  to  obtain,  since  these  equations  are  too  complex 
and  have  high  coupling.  In  addition,  Eq.  (4.30)  has  a  singular  point  at  rj  = 
0.  Moreover,  as  pqinted  out  by  Singleton  [17],  regardless  of  the  choice  of 
u'  '  (0)  and  (0),  the  resulting  solutions  to  u^'(n)  and  T^'(r))  from 

Eqs.  (4.28)  and  (4.29)  always  approach  zero  as  n  approaches  infinity,  making 
it  impossible  to  pick  out  the  correct  solutions.  Therefore,  in  this 
analysis,  it  was  not  attempted  to  obtain  the  first-order  solutions  to  Eqs. 
(4.27)-(4.30)  but  just  the  zeroth-order  solutions  to  Eqs.  (4.21)-(4. 24 ) 
which  are  of  more  practical  interest.  The  first-order  problem  for  the 
particle  velocity  and  temperature  are  readily  obtained.  With  the  series 
expansion  (4.20),  the  slip  quantities  of  first  order  are  given  as 


(1) 


_  f 


(o)u(or 


(4.32) 


2T 


(0) 


w 


(1) 


1 


(nu«»2  ♦  -  li£ll  f(°)2)  (4.33) 

2j(0) w  T<°> 


.(1) 


_  3Pr  f^ua 


(0)T(0)' 


4a 


-(0) 


(4.34) 


Equat i ons 
values 


is  ( 4 . 32 ) - ( 4 . 34 )  are  not  differential  but  alaebraic  equati 
of  three  slip  quantities  u^',  v^  and  at  any  given 


s  *  s  s 

determined  from  the  zeroth-order  solutions  for  the  gas  phase. 


ons.  The 
point  are 


5.0  RESULTS  AND  DISCUSSIONS 


The  zeroth  and  first-order  equations  for  the  large-slip  and  small-slip 
limits  can  be  solved  numerically.  They  are  a  system  of  nonlinear, 
second-order ,  ord i nary-di f ferent i al  equations.  The  corresponding  boundary 
conditions  are  specified  at  the  two  end  points,  i.e.,  at  the  wall  and  the 
outer  edge  of  the  boundary  layer.  Mathematically,  it  is  a  two-point 
boundary-val ue  problem  and  it  can  be  solved  by  Gear's  method  [22]. 
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From  Eqs.  (3. 18 ) - ( 3 . 20 ) ,  the  zeroth-order  problem  for  the  gas  phase  in 
the  large-slip  limit  is  as  simple  as  that  for  the  boundary-layer  flow  of  a 
pure  gas  without  particles.  Therefore,  as  for  the  conventional  viscous 
flows  of  a  pure  gas,  the  Reynolds  number,  the  Prandtl  number,  the  Eckert 
number  and  the  viscosity  power  index  are  important  controlling  parameters  in 
the  analysis  of  compressible,  laminar,  boundary-1 ayer  flows  of  a 
gas-particle  mixture.  The  numerical  solutions  of  Eqs.  (3.18)-(3.20)  are 
given  in  Figs.  3  to  5  and  the  influence  of  the  parameters  Pr,  Ec  and  won 

the  flow  properties  are  shown  in  these  figures  and  it  is  seen  that  their 

effects  are  relatively  small.  Here,  it  is  not  necessary  to  discuss  the 
zeroth-order  solutions  for  the  gas  phase  in  detail  since  it  is  the  same  as 
the  similarity  solution  for  the  flat-plate  boundary  layer  of  a  pure  gas. 
Similarly,  the  zeroth-order  solutions  for  the  particle  phase  are  readily 
obtained.  Equations  ( 3 .21 ) -( 3 .23)  indicate  that,  in  the  zeroth-order 
approximation,  the  particle  motion  in  the  boundary  layer  remains  uniform. 
All  the  zeroth-order  flow  quantities  for  the  particle  phase  (density, 
velocity  and  temperature)  are  the  same  as  those  in  the  freestream  or  the 

external  flow.  This  is  due  to  the  fact  that  both  the  gas  and  the  particles 

move  independently  of  each  other  in  the  zeroth-order  problem.  The  influence 
of  the  particles  on  the  flow  properties  is  prevalent  only  in  the  first  or 
higher  order  solutions.  It  is  a  major  feature  of  the  two-phase 
boundary-layer  flows  in  the  large-slip  region. 


Figures  6  to  8  show  the  first-order  solutions  for  the  gas  phase,  i.e., 
the  solutions  to  Eqs.  (3.37)-(3. 39 ) .  They  are  the  numerical  results  for  the 
Stokes  case  and  the  effects  of  Ec,  Pr  and  w  are  significant  this  time.  For 
the  first-order  problem,  the  same  value  of  the  flow  parameters  Pr,  Ec ,  wand 
Tw  were  chosen  as  in  the  zeroth-order  problem  where  Pr  =  0.69-1.0,  Ec  = 
0. 1-1.0,  w  =  0. 5-1.0  and  Tw  =  0.5.  For  the  non-Stokes  case,  the  numerical 
results  with  Pr  =  0.69,  Ec  =  1.0,  w  =  0.67  and  Tw  =  0.5  are  presented  in 
Fig.  9.  They  cover  quite  a  wide  range  of  the  particle  Reynolds  number 
(Rep  =  0.1-100.0)  and  the  changes  are  very  significant.  For  the  particle 
phase,  the  first-order  solution  can  be  obtained  by  numerically  integrating 
Eqs.  ( 3 .46 ) - ( 3 .49 ) .  The  results  for  the  Stokes  case  with  a  =  1.0  are  given 
in  Figs.  10  to  12.  It  is  seen  that;  significant  changes  occur  in  pp  and  T^1 ' 
with  w,  Pr  and  Ec  and  in  fp  with  Pr  and  Ec.  The  results  for  the 
non-Stokes  case  are  shown  in  Fig.  13  where  the  changes  with  Rep  are  even 
more  significant.  The  computations  for  the  particle  phase  is  carried  out 
under  the  same  conditions  as  those  for  the  gas  phase. 


By  comparing  the  results  for  the  non-Stokes  case  with  those  for  the 
Stokes  case,  it  is  seen  that  the  results  based  on  the  Stokes  relation  are 
reasonable  qual i tat i vely .  They  present  the  similar  tendency  of  variations 
in  the  flow  properties,  such  as  velocity  and  temperature.  However,  they  are 
not  correct  quant i tat i vely ,  especially  for  the  large  particle  Reynolds 
number,  as  expected.  Nevertheless,  this  comparatively  simple  case  of  the 
Stokes  relation  is  still  considered  in  many  analyses  of  dusty-gas  flows 


since  it  is  useful  for  understanding  the  main  characteristics  of  two-phase 
flow  phenomena. 

From  the  solution  to  the  first-order  problem,  it  was  found  that  there 
exist  significant  differences  in  the  first-order  flow  profiles  between  the 
two  phases.  For  instance,  the  first-order  velocity  of  the  gas  is  positive 
and,  while  passing  across  the  boundary  layer  from  the  outer  edge  to  the 
wall,  it  increases  first  to  a  maximum  value  and  then  decreases  to  zero  (see 
Fig.  7).  By  rontrast,  the  first-order  velocity  of  the  particles  is  negative 
and  its  magnitude  increases  monotonical ly  from  the  outer  edge  of  the 
boundary  layer  to  the  wall  (c.f.  Fig.  11).  This  arises  from  the  fact  that 
the  mechanisms  of  motion  for  the  two  phases  are  not  the  same.  There  are  two 
kinds  of  forces  exerted  on  the  gas:  the  viscous  force  by  the  gas  and  the 
drag  force  by  the  particles.  For  the  particle  phase,  however,  only  the  drag 
force  of  the  gas  influences  its  motion.  Therefore,  after  entering  the 
boundary  layer  at  the  leading  edge,  the  gas  decreases  immediately  its 
tangential  velocity  from  the  freestream  value  at  the  outer  edge  to  zero 
velocity  at  the  wall  due  to  viscosity.  Since  the  density  of  the  particle 
material  is  much  greater  than  the  gas  density,  the  particles  cannot 
accommodate  this  rapid  deceleration  but  tend  to  slip  through  the  gas  as  they 
decelerate.  It  takes  some  time  for  the  particles  and  gas  to  adjust  to  an 
equilibrium  state.  It  implies  that  in  the  large-slip  region  near  the 
leading  edge,  the  gas  has  small  deviation  from  the  pure-gas  boundary-layer 
flow  while  the  particles  have  small  deviation  from  their  original  state  of 
uniform  motion  in  the  freestream.  The  particles  are  'frozen'.  This 
situation  is  justified  by  the  zeroth-order  solution,  which  represents  the 
complete  frozen-flow  limit.  The  relaxation  process  takes  place  throughout 
the  equilibrium  length  \*.  In  the  meantime,  owing  to  the  slip  velocity,  the 
drag  force  arises  between  the  two  phases  and  then  the  first-order  flow  is 
induced  by  this  gas-particle  interaction.  The  gas  is  accelerated  and  the 
particles  are  retarded.  This  is  the  reason  why  the  two  phases  have  their 
first-order  velocities  in  opposite  directions.  While  traversing  the 
boundary  layer  from  the  outer  edge  to  the  wall,  the  slip  velocity  increases 
and  then  the  first-order  velocity  for  the  two  phases  both  increase  first  in 
the  region  near  the  outer  edge,  since  the  drag  force  is  proportional  to  the 
slip  velocity.  The  first-order  velocity  of  the  particles  continues  to 
increase  in  magnitude  on  approaching  the  wall,  since  the  particle  motion  is 
driven  only  by  the  drag  force.  However,  for  the  gas  phase,  in  the  region 
near  the  wall  where  the  velocity  gradient  for  the  gas  is  great,  the  viscous 
force  prevails  and  the  no-slip  condition  at  the  wall  forces  the  gas  velocity 
to  go  to  zero.  Thus,  the  first-order  velocity  of  the  gas  decreases  in  the 
inner  boundary  layer  and  vanishes  at  the  wall.  A  similar  argument  is  valid 
for  the  first-order  temperature  profit’  by  employing  the  correspondence  of 
the  temperature  to  the  velocity,  the  heat  conductivity  to  the  viscosity  and 
the  heat  transfer  to  the  drag  force. 

Finally,  for  boundary- 1 ayer  analyses,  there  are  three  characteristic 
quantities  of  interest:  the  shear  stress  at  the  wall  the  heat-transfer 
rate  at  the  wall  q*  and  the  displacmeent  thickness  &*.  As  usual,  they  are 


determined  from  the  flow  profiles  of  the  gas  phase: 
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It  is  convenient  to  introduce  the  following  nondimensional  characteristic 
quantities: 
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where  Re^  is  the  flow  Reynolds  number  based  on  the  freestream  velocity  u* 
and  the  velocity-equilibrium  length  X.*, 
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Then,  the  nondimensional  boundary-1 ayer  characteri sti cs  can  be  expressed 
as: 
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where  is  the  nondimensional  viscosity  of  the  gas  at  the  wall 


and  the  coefficients  aj,  a2,  blf  b2,  Cj  and  c2  are  given  by 


a  l  =  u<°)' (0), 


bx  =  t(°)‘(0). 


a  =  u 

a2 


1 ^ 1 (0) 
u(°)‘(0) 


(5.8) 


=  T(1)l (0) 
t(°) ' (0) 


(5.9) 


c‘  *  J0  (1'7ior)d" 


u(0)T(i)  -  TWam  d. 


(5.10) 


From  the  above  relations,  it  is  known  that  the  coefficients  a j ,  bj  and  Cj 
are  determined  only  by  the  zeroth-order  solution  and  that  the  coefficient 
a2,  b2  and  c2  depend  on  the  first-order  solution  as  well  as  the  zeroth-order 
solution.  In  fact,  the  first  three  coefficients,  i.e.,  ax,  bt  and  cL,  give 
the  zeroth-order  approximation  of  the  three  characteristics  which  is  the 
same  as  for  the  similarity  solution  of  a  pure  gas.  The  other  three 
coefficients  a2,  b2  and  c2  represent  the  first-order  modification  owing  to 
the  presence  of  the  particles.  These  coefficients  can  be  estimated  from  the 
numerical  results  using  eqs.  (5.8)-(5. 10) .  In  Table  1,  the  listed  values 
are  the  coefficients  for  the  case  where  the  flow  parameters  are  Pr  =  0.69, 
Ec  =  1.0,  u  =  0.67  and  Tw  =  0.5. 


Table  1 

Coefficient  Values  for  Boundary-Layer  Characteristic  Quantities 


Stokes' 

non-Stokes'  Case 

La  Se 

Re=0. 1 

Re=1.0 

Re=10.0 

Re=100.0 

0.5472 

0.5472 

0.5472 

0.5472 

0.5472 

1.888 

1.599 

2.312 

3.728 

9.870 

0.4382 

0.4382 

0.4382 

0.4382 

0.4382 

2.275 

2.058 

2.827 

4.487 

11.21 

1.101 

1.101 

1.101 

1.101 

1.101 

-0.6672 

-0.4515 

-0.7152 

-1.148 

-2.748 

From  the  values  given  in  the  above  table,  it  is  seen  that  the  shear  stress 
and  heat-transfer  rate  at  the  wall  in  the  case  of  dusty  gases  become  greater 
than  those  in  the  case  of  pure  gases  and  the  displacement  thickness 
thi nner. 


5.2  Small-Slip  Limit 

The  zeroth-order  equations  for  the  gas  as  well  as  for  the  particles  in 
the  small-slip  limit,  Eqs.  (4.21)-(4.23) ,  are  similar  to  the  conservation 
equations  for  a  pure-gas  boundary-layer  flow  but  with  modified  properties. 
Physically,  the  small-slip  approximation  represents  a  quasi -equi 1 ibri urn  flow 
and  the  zeroth-order  problem  constitutes  the  exact  equilibrium  limit  where 
the  particles  are  'fixed'  to  the  mass  of  the  gas  so  that  the  gas  and 
particles  move  together  like  a  perfect-gas  mixture.  For  the  dilute 
two-phase  system  with  the  approximation  of  negligible  volume  fraction  of  the 
particles,  the  particles  contribute  to  the  mixture  density  but  not  to  the 
viscosity  [23]: 


(1+P)p* 


(5.11) 


p*  =  p* 


(5.12) 


The  other  thermodynamic  properties  are  given  by 


„*  _  1  +  P/a 

P  CP  1  +  p 


(5.13) 


k*  =  k* 


(5.14) 


Then  the  modified  similarity  parameters  can  be  expressed  as 


tS1 
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c*  T* 
p  00 
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Substituting  the  modified  Prandtl  number  and  Eckert  number  into  Eqs. 
(3.18)-(3.20) ,  which  are  exactly  the  same  as  the  boundary-layer  equations 
for  pure  gases  as  mentioned  before,  the  resulting  'modified'  equations  are 
just  the  zeroth-order  equations  for  the  small-slip  limit,  Eqs. 
(4.21)-(4.23) .  It  implies  that  in  the  zeroth-order  approximation,  the 
gas-particle  mixture  behaves  like  a  pure  gas  with  modified  thermodynamic 

properties.  In  this  paper,  the  numerical  solutions  in  the  small-slip  limit 
are  calculated  for  the  condition  of  a  =  1.0.  Under  this  condition,  the 
zeroth-order  equations  (4.21)-(4.23)  for  the  small-slip  limit  reduce  to 

those  for  the  large-slip  limit,  Eqs.  ( 3 . 18) - ( 3 . 20 ) .  The  results  are  given 
in  Figs.  3-5.  Clearly,  in  the  small-slip  limit,  the  zeroth-order  velocity 
and  temperature  for  the  particles  are  the  same  as  the  ones  for  the  gas. 
From  Eq.  (4.24),  it  is  found  that  the  nondimensional  density  of  the  particle 
phase  p(°'  is  equal  to  the  nondimensional  density  of  the  gas  phase  p(°'  (or 
the  reciprocal  of  the  nondimensional  gas  temperature  1/t(°').  It  is  seen 
from  Fig.  14  that  p(^',  or  p'^',  varies  monotonical ly  from  its  maximum  value 
at  the  wall  to  its  freestream  value  at  the  outer  edge.  In  this  report,  the 
constant  wall  temperature  Tw  is  specified  as  Tw  =  0.5  and  then  the 

density  at  the  wall  is  equal  to  2,  as  shown  in  Fig.  14.  From  Eqs.  (4.9)  and 

(4.10),  the  densities  in  the  dimensional  form  are  given  by 


P*  =  P®  P,  Pp  =  Pp®  Pp  (5.17) 

Neglecting  the  small  quantities  of  first  order,  the  above  gas  and  particle 
densities  are  approximated  by 

P*  =  PS,  p(°),  Pp  =  Pp®  Pp  ^  (5.18) 

Therefore,  at  all  points  of  the  boundary  layer  in  the  small-slip  region, 

i£  =  p  (5.19) 

P* 

It  means  that  the  constant  loading  ratio  of  the  particles  holds  across  the 
whole  boundary  layer  in  the  small-slip  region.  In  other  words,  the  solid 
particles  remain  attached  to  their  original  gas  mass  and  always  move 
together  with  this  gas  mass.  The  two-phase  system  behaves  like  a  gaseous 
mixture.  It  is  a  major  feature  of  the  gas-particle  flow  in  the  small-slip 
region. 

The  slip  quantities  u^1^,  v^  and  T^1^  are  given  in  Figs.  15  to  17, 
where  the  effects  of  the  flow  parameters  Pr,  Ec  and  w  on  the  first-order 


flow  of  the  particles  are  shown.  In  fact,  the  slip  quantities  represent  the 
first-order  approximation  for  the  particle  phase.  It  is  seen  that  the 
profile  of  the  normal  slip  velocity  is  different  from  that  of  the 

tangential  slip  velocity  u^'.  At  the  outer  edge  of  the  boundary  layer,  the 
tangential  slip  velocity  becomes  zero  but  the  normal  slip  velocity 
approaches  a  finite  value,  since  the  boundary  conditions  for  the  tangential 
and  normal  velocities  are  different  in  boundary- 1 ayer  analyses.  As  in  the 
usual  boundary-1 ayer  problems,  the  tangential  particle  velocity  at  the  outer 
edge  of  the  boundary  layer  should  be  equal  to  the  freestream  value  and  then 
the  slip  velocity  should  become  zero  at  the  outer  edge.  However,  no  similar 
boundary  conditions  at  the  outer  edge  can  be  specified  for  the  normal 
velocity.  The  unique  boundary  condition  for  the  normal  velocity  is  that  it 
is  equal  to  zero  at  the  wall.  Owing  to  the  continuity  equation,  the  normal 
velocity  is  induced  and  approaches  its  maximum  value  at  the  outer  edge. 
Therefore,  the  normal  slip  velocity  at  the  outer  edge  takes  a  finite  value 
which  is  the  difference  between  the  normal  velocities  of  the  two  phases  at 
the  outer  edge.  In  addition,  by  comparing  Fig.  15  with  Fig.  11,  it  is  found 
that  there  exist  significant  changes  in  the  first-order  profile  for  the 
tangential  velocity  of  the  particles  in  the  two  limiting  regions.  In  the 
small-slip  region,  the  particle  slip  velocity  at  the  wall  is  equal  to  zero, 
while  in  the  large-slip  region,  the  first-order  velocity  of  the  particles 
has  its  maximum  value  in  magnitude  at  the  wall  as  mentioned  earlier.  As  a 
result  of  the  maximum  slip  velocity,  the  interaction  term  between  the  gas 
and  the  particles  has  its  maximum  value  at  the  wall  and  then  the  maximum 
deceleration  of  the  particles  takes  place  along  the  wall.  Hence,  at  some 
distance  from  the  leading  edge,  the  particle  velocity  at  the  wall  reduces  to 
zero  and  is  equal  to  that  of  the  gas.  After  this  point  the  particles  keep 
their  zero  velocity  at  the  wall  because  of  the  zero  slip  velocity.  This 
special  point,  where  the  particle  velocity  becomes  zero  at  the  wall,  is 
defined  as  the  critical  point  for  the  gas-particle  boundary-layer  flow.  At 
the  critical  point,  the  dusty-gas  boundary  layer  essentially  fulfills  the 
transition  from  the  quasi-frozen  flow  to  the  quasi -equi 1 i bri urn  flow.  The 
two-phase  flow  in  the  small-slip  limit  is  a  typical  example  of  a 
quasi -equi  1  i  bri  urn  flow.  By  comparing  Fig.  17  with  Fig.  12,  the  same 
situation  happens  to  the  first-order  temperature  of  the  particles  in  the  two 
limiting  regions:  at  the  wall,  the  first-order  temperature  has  its  maximum 
value  in  the  large-slip  region  and  the  temperature  defect  vanishes  in  the 
small-slip  region.  A  detailed  discussion  is  omitted  here,  since  it  is 
similar  to  the  above  case  of  velocity.  Similar  to  the  large-slip  case,  it 
is  interesting  to  obtain  the  expressions  for  the  bounda ry- 1  aye r 
characteristic  quantities.  However,  in  the  small-slip  case,  since  there  are 
no  solutions  available  to  the  first-order  equations  of  the  gas,  they  can  be 
expressed  only  in  zeroth  order.  Similarly,  the  three  boundary-1 ayer 
characteristics  are  given  in  the  nondimensi onal  form: 


) 

i 
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qw  =  -  —  M  +  P  bj)  (5.21) 

✓x 


6  =  /x  (4L)  (5.22) 
/1+p 

The  coefficients  a^,  bj  and  Cj  in  Eqs.  (5.20)-(5.22)  have  the  same  values  as 
in  Table  1.  By  contrast  with  the  large-slip  limit,  the  zeroth-order 
expressions  for  the  boundary-layer  characteri sties  in  the  small-slip  limit, 
Eqs.  ( 5 . 20) -(5.22) ,  involve  the  effects  of  the  particles.  In  fact,  the  term 

/1+P  in  Eqs.  ( 5 . 20 ) - (5 .22)  represents  the  alteration  of  the  boundary  layer 
by  the  particles.  In  the  small-slip  limit,  the  two-phase  system  acts  like  a 
single  gaseous  system  with  modified  properties  as  pointed  out  before.  With 
Eqs.  (5.9)-( 5.10) ,  the  'modified'  similarity  variable  r)  becomes 


f~ :i+p)u* 

/  2  v*“ x* 


j*  = 


h 


/1+p 


(5.23) 


This  implies  that  the  boundary-layer  flow  of  a  dusty  gas  in  the  small -si ip 
limit  corresponds  to  a  similarity  solution  with  the  normal  scale  modified  by 

the  factor  /1+p,  owing  to  the  particles.  Consequently,  the  shear  stress  and 
heat-transfer  rate  at  the  wall  increase  and  the  displacement  thickness 

decreases  by  the  same  factor  of  /1+  p.  Therefore,  it  can  be  concluded  that 
the  presence  of  particles  enhances  the  shear  stress  and  heat-transfer  at  the 
wall  and  thins  the  boundary  layer  in  the  two  limiting  regions.  This 
tendency  can  be  seen  in  Figs.  18  and  19,  where  the  shear  stress  and 

heat-transfer  rate  are  shown  as  functions  of  the  nondimensional  distance  x 

for  the  cases  with  and  without  particles.  As  expected,  the  results  for  the 

large-slip  limit  and  the  small-slip  limit  coincide  with  the  pure-gas  results 
in  the  limits  x  +  0  and  x  +  ®,  respectively.  Note  that  the  large-slip 
results  when  x  >  0.1  and  the  small-slip  results  when  x  <  10  are  meaningless 
since  the  aysmptotic  solutions  are  not  valid.  Physically,  the  changes  in 
the  characteristics  caused  by  the  particles  can  be  explained  as  follows. 
The  gas-flow  profiles  with  and  without  particles  are  schematically  shown  in 
Fig.  20.  As  a  result  of  the  interaction,  the  gas  velocity  and  temperature 
increase  in  the  cold-wall  case  (say,  Tw  =  0.5).  Then  the  derivatives  of 
the  gas  velocity  and  temperature  with  respect  to  the  normal  coordinate  y*  at 
the  wall  become  greater  than  those  without  particles.  These  changes  result 
in  an  increase  in  the  shear  stress  and  heat-transfer  at  the  wall,  since  they 
are  proportional  to  those  derivatives.  In  addition,  the  boundary  layer 
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becomes  thinner  since  the  velocity  approaches  its  freestream  value  more 
qui ckly . 

The  numerical  results  for  the  asymptotic  solutions  using  the 
series-expansion  method  were  compared  with  the  finite-difference  solutions 
in  the  two  limiting  regions.  The  agreement  between  the  asymptotic  solutions 
and  the  difference  solutions  was  excellent  [1]. 


6.0  CONCLUDING  REMARKS 


Some  general  conclusions  obtained  from  the  asymptotic  solutions  of  the 
flat-plate  boundary-layer  flow  of  a  dilute  gas-particle  mixture  are 
summarized  as  follows: 


(1)  The  asymptotic  solutions  to  the  dusty-gas  boundary-layer  equations  can 
be  obtained  using  a  series-expansion  method.  They  describe  the 
limiting  properties  of  two-phase  flows  in  the  large-slip  region  and  the 
small-slip  region,  which  are  characteri zed  by  a  frozen  flow  and  an 
equilibrium  flow,  respectively.  The  asymptotic  solutions  are  in 
excellent  agreement  with  the  finite-difference  solutions. 

(2)  The  interaction  between  the  gas  and  particles  determines  the  flow 
properties  of  the  particle  phase,  and  influences  strongly  the  flow 
properties  of  the  gas  phase  in  addition  to  the  viscosity.  When  the 
particle  slip  Reynolds  number  is  high,  a  proper  expression  for  the  drag 
and  heat  transfer  between  the  two  phases  should  be  specified  instead  of 
using  Stokes'  relations.  The  results  when  using  Stokes'  relations  are 
reasonable  qualitatively  but  not  correct  quantitatively  for  a  Reynolds 
number  greater  than  unity. 

(3)  For  a  given  gas-particle  system  with  specified  values  of  the  mass 
loading  ratio  and  the  ratio  of  the  specific  heats  of  the  two  phases, 
similar  to  the  case  of  high-velocity  viscous  flows  of  a  pure  gas,  all 
of  the  following  parameters  are  important  for  the  analysis  of 
compressible  laminar  boundary-layer  flows  of  gas-particle  mixtures: 
Reynolds,  Prandtl  and  Eckert  numbers  as  well  as  the  transport 
properties  (viscosity  and  heat  conductivity). 


(4) 


For  compressible,  laminar,  boundary-layer  flows  of  dusty  gases,  the 
shear  stress  and  the  heat-transfer  rate  at  the  wall  increase  and  the 
displacement  thickness  decreases  when  compared  with  the  corresponding 
results  for  a  pure  gas.  Owing  to  the  presence  of  the  particles,  the 
gas  velocity  and  temperature  increase  on  a  flat-plate  boundary  layer 
with  a  cold  wall.  As  a  result,  the  velocity  and  temperature  gradients 
at  the  wall  for  the  gas  phase  increase  so  that  the  shear  stress  and 
heat  transfer  are  enhanced  and  the  velocity  achieves  its  freestream 
value  at  a  shorter  distance  from  the  wall  so  that  the  displacement 
thickness  of  the  boundary  layer  is  decreased. 
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APPENDIX 


RELAXATION  PROCESS  AND  VELOCITY  EQUILIBRIUM  TIME 


In  gas-particle  flows,  generally  speaking,  the  gas  and  the  particles 
may  have  different  velocities  and  temperatures,  and  as  an  immediate 
consequence  the  two  phases  must  interact.  Every  particle  experiences  a  drag 
force  and  exchanges  heat  with  the  gas.  The  same  happens  to  the  gas  but  in 
the  opposite  direction.  Consequently,  the  velocities  and  temperatures  tend 
to  approach  each  other,  and  the  instantaneous  rate  of  this  approach  depends 
on  the  instantaneous  velocity  and  temperature  differences  between  the 
particle  and  the  gas.  This  phenomenon,  which  is  usually  known  as  the 
relaxation  process,  is  a  major  feature  of  gas-particle  flows.  A  discussion 
of  the  viscous  drag  between  the  gas  and  the  particles  and  of  the  velocity 
relaxation  process  will  illustrate  the  relaxation  processes  in  a  two-phase 
system. 

For  the  case  of  a  single  spherical  particle  in  a  viscous  flow,  if  the 
particle  velocity  Up  (here  the  unstarred  quantities  denote  dimensional  flow 

properties)  is  different  from  the  gas  velocity  U,  the  viscous  drag  exerted 

by  the  gas  on  the  particle  depends  on  the  relative,  or  slip,  velocity 
•>  •>  -> 

(Up  -  U).  As  usual  in  fluid  mechanics,  the  drag  F  is  expressed  by  means  of 
a  drag  coefficient  Cy  which  is  defined  in  terms  of  the  dynamic  head  of  the 
relative  flow  and  the  frontal  area  of  the  particle: 


F  =  CD{I  p(U  -  Up)  | U  -  Up |  |d2}  ( A.  1 ) 

where  p  is  the  gas  density  and  d  is  the  particle  diameter.  On  the 
right-hand  side  of  Eq.  (A.l),  the  square  of  the  relative  velocity  is  written 
in  this  manner  to  ensure  that  the  drag  force  always  has  the  correct  sign. 
The  drag  coefficient  Cq  is  a  function  of  the  particle  slip  Reynolds 
number 

CD  =  C0(Res)  (A. 2) 

where  Res  is  the  slip  Reynolds  number  based  on  the  particle  diameter 


Re 


s 


(A. 3) 


A.l 


where  ii  is  the  gas  viscosity. 


For  the  Stokes  case,  with  a  low-slip  Reynolds  number  less  than  about 
one,  the  drag  coefficient  takes  the  simple  form 


(A. 4) 


This  is  the  Stokes  relation  for  the  drag  force.  For  a  higher  slip  Reynolds 
number  (Res  >  1),  there  are  several  empirical  relations  for  the  drag 

coefficient  given  by  different  investigators.  It  is  convenient  to  write 


CD  =  CD0  *  D(Res> 


D  = 


(A. 5) 


(A. 6) 


where  D  is  the  so-called  normalized  drag  coefficient  in  this  analysis, 
is  clear  that  D(Res)  =  1  corresponds  to  Stokes'  drag. 


The  above  discussion  deals  only  with  isolated  particles.  Of  course,  if 
there  is  an  appreciable  number  density  of  particles,  the  effective,  or 
apparent,  drag  coefficient  may  be  different  from  that  for  a  single  particle 
because  of  some  forms  of  interaction  between  the  particles,  such  as  direct 
collisions  and  particle-wake  interaction.  At  present,  there  are  no  adequate 
analytical  or  experimental  results  available  for  these  interactions.  In 
this  report,  the  case  of  a  dilute  gas-particle  mixture  is  considered,  and  it 
is  reasonable  to  assume  that  the  drag  coefficient  for  a  single  sphere  is 
valid  for  the  particle  cloud.  Therefore,  if  the  drag  coefficient  takes  the 
form  of  the  Stokes  relation,  the  drag  force  on  every  particle  is 


v^-v|d2} 


24  l*d2  _  iT  m.T 


{S_  p(U  -  UD)|U  -  U_|  } 


p(Up-U|d/n 


3KMd(U  -  Up) 


(A. 7) 


>  *-  JV  '»■ 


'  VA'.V’.'A  1 


VWV' 


Let  rip  denote  the  number  density  of  the  particles,  pp  the  density 
of  the  particle  phase,  ps  the  density  of  the  particle  material  and  m  the 
mass  of  a  particle.  Then 


(A. 8) 


n„  =  = 


(A. 9) 


The  force  per  unit  volume  of  the  gas-particle  mixture  acting  on  the  gas, 
■>  — 

D  =  iDx  +  jDy,  can  be  given  as 


D  =  -np  •  F 


- -  •  3npd(Un  -  U) 

o  •  11  d3 
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U„  -  U 


(A. 10) 


where  tv  is  defined  as  the  particle  velocity-equilibrium  time 


(A. 11) 


The  velocity-equilibrium  time  tv  depends  on  the  density  of  particle 
material,  the  particle  diameter  and  the  gas  viscosity.  Some  typical  values 
of  the  velocity-equilibrium  time  are  given  in  Table  2,  where  the  particle 
material  density  pp  =  2.5  g/cm3  (typical  of  crown  glass). 


Table  2 


Typical  Values  of  the  Particle  Velocity  Equilibrium  Time  (ms) 


Gas 

Temperature 

T  (K) 

Particle  Diameter  d 

( pm) 

10 

20 

40 

300 

0.755 

3.02 

12.08 

400 

0.605 

2.42 

9.68 

500 

0.510 

2.04 

8.16 

Physically,  xv  is  the  time  required  for  a  particle  to  reduce  the  slip 
velocity  to  e-1  of  its  initial  value.  To  examine  this  criteria,  consider  a 
flow  field  with  constant  gas  velocity  and  temperature.  For  the  particles 
moving  in  this  flow  field,  Newton's  law  of  motion  becomes 


m 


dt 


(A. 12) 


Substituting  Eqs.  (A. 7)  and  (A. 8),  Eq.  (A. 12)  becomes 


5>  .i 


dt 
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3upd(U  -  Up) 
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=  lit-  (U  -  U  ) 
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(A. 13) 


A. 4 


Since  the  gas  velocity  and  temperature  are  assumed  as  constant,  i.e.,  U  and 
ty  are  constant,  Eq.  (A. 13)  is  readily  integrated  with  the  result 


-  U 


(U 


p,o 


(A. 14) 


where  Up  Q  is  the  initial  velocity  of  the  particles  at  t  =  0.  This  result 
indicates’  a  typical  velocity-relaxation  process  in  which  the  slip  velocity 
decays  exponentially.  From  Eq.  (A. 14),  it  is  known  that  the  velocity- 
equilibrium  time  tv  is  the  time  elapsed  for  a  particle  to  reduce  its 
original  relative  velocity  by  e-1.  It  should  be  pointed  out  that,  in  fact, 
the  velocity  relaxation  process  is  not  completed  in  the  velocity-equilibrium 
time  tv ,  since  the  slip  velocity  still  has  an  appreciable  value,  i.e., 
about  37%  of  its  original  value.  In  addition,  if  the  gas  velocity  and 
temperature  are  not  constant,  or  Stokes'  drag  does  not  apply,  Eq.  (A. 13) 
becomes  more  difficult  to  solve,  and  the  solution  is  no  longer  a  simple 
exponential  decay  of  the  slip  velocity,  Eq.  (A. 14).  Nevertheless,  the 
parameter  tv  remains  a  convenient  measure  for  the  tendency  of  the 
particles  to  reach  velocity  equilibrium  with  the  gas. 

For  the  non-Stokes  case,  the  drag  force  on  every  particle  is  given  by 


F  •  CD()|i  P(U  -  ffpHU-  Up|  |d2)D 


=  3rcud(U  -  Up)D 


(A. 15) 


The  drag  force  per  unit  volume  exerting  on  the  gas  by  the  particles  is 
obtained  as 
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where  u  and  v  are  the  velocity  components  for  the  gas,  and  uD  and  vp  for 
the  particles: 

U  =  iu  +  jv,  Up  =  iup  +  j  vp  (A. 18) 

The  dissipation  function  due  to  the  relative  motion  of  the  particles  is 

♦  =  D  -  (Up  -  U) 

=  (up  -  u)Dx  +  (vp  -  v)Dy  (A. 19) 

The  heat  transfer  from  a  particle  to  the  gas  has  the  form 


where  k  is  the  heat  conductivity  of  the  gas,  Nu  is  the  Nusselt  number,  and 
Tp  and  T  are  the  particle  and  gas  temperatures,  respectively.  Then  the 
total  heat  transfer  to  the  gas  per  unit  volume  is  given  by 


(A. 20) 


where  c$  is  the  specific  heat  of  the  particle  phase  and  tj  is  the 
temperature-equilibrium  time: 


The  temperature-equilibrium  time  xy  can  be  expressed  in  terms  of  the 
velocity-equilibrium  time  xv: 


It  =  Ps<j2  18  P 

T  18  p  12k 


c 


s 


(A. 22) 


where  Cp  is  the  specific  heat  at  constant  pressure  for  the  gas  and  Pr  is 
the  gas  Prandtl  number: 


Pr 


-fRlt 
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(A. 23) 


FIGURE  1.  SCHEMATIC  OF  A  DUSTY-GAS  BOUNDARY-LAYER  OVER  A  SEMI-INFINITE  FLAT 
PLATE.  REGION  I:  LARGE-SLIP  REGION.  REGION  II:  MODERATE-SLIP 
REGION.  REGION  III:  SMALL-SLIP  REGION. 
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FIGURE  3  -  CONTINUED 
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ZEROTH-ORDER  VELOCITY  u<°)  FOR  THE  GAS  PHASE. 

Pr  =  1.0,  Ec  =  1.0  AND  w  =  0.5,  0.67,  0.  77  AND  1  0 

RESPECTIVELY.  '  ‘ 

Ec  =  1.0,  tu  =  0.67  AND  Pr  =  0.69,  0.  72,  0.  75  AND  1  0 

RESPECTIVELY.  ’  ’ 

Pr  =  0.69,  io  »  0.67  AND  Ec  =  0.1,  0.4  AND  1.0,  RESPECT  I  VtLY. 
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FIGURE  5.  THE  ZEROTH-ORDER  TEMEPRATURE  T^0)  FOR  THE  GAS  PHASE. 

(a)  Pr  =  1.0,  Ec  =  1.0  AND  u  =  0.5,  0.67,  0.  77  AND  1.0 

RESPECTIVELY. 

(b)  Ec  =  1.0,  w  =  0.67  AND  Pr  =  0.69,  0.  72,  0.  75  AND  1.0 

RESPECTIVELY. 

(c)  Pr  =  0.69,  oo  =  0.67  AND  Ec  =  0.1,  0.4  AND  1.0,  RESPECTIVELY. 
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FIGURE  8  -  CONTINUED 
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FIGURE  9.  FIRST-ORDER  SOLUTIONS  FOR  THE  GAS  PHASE  IN  THE  LARGE-SUP  LIMIT 
FOR  THE  NON-STOKES  CASE  WITH  Pr  =  0.69,  Ec  =  0.67  ANU  Re  =  0.1, 
1.0,  10.0  AND  100.0,  RESPECTIVELY.  V 

(a)  FUNCTION  F(1);  (b)  VELOCITY  U( 1 >  (c)  TEMPERATURE  ®( 1 > . 


vf!< 

,V#5 


m 


vvl 

•K 


$ 


0'"  .  H 

1  .  0 


1. RE-O. 100 

2.  RE-1. 000 

3.  RE-10. 00 

4.  RE-100.0 


co-0. 666 
PR- 0.690 
EC-1.000 


p  n  L  _  _  i_  _  -  -L  —  —  -1  --1 - 1 - 1 - 1 - « - 1 

O.c  l.n  3.0  4.0  6.0  6.0  7.0  8.0  8.0  ]0.n 


FIGURE  9  -  CONTINUED 
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FIGURE  10.  FIRST-ORDER  DENSITY  pt1)  FOR  THE  PARTICLE  PHASE  IN  THE 


LARGE-SLIP  LIMIT  FOR  THE  STOKES  CASE. 
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FIGURE  11  -  CONTINUED 
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FIGURE  12.  FIRST-ORDER  TEMPERATURE  T^1)  FOR  THE  PARTICLE  PHASE  IN  THE 
LARGE-SLIP  LIMIT  FOR  THE  STOKES  CASE. 

(a)  Pr  =  1.0,  Ec  =  l.n  AND  uj  -  0.5,  0.67,  0.  77  AND  1  0 
RESPECTIVELY. 

(b)  Ec  -  1.0,  w  -  0.67  AND  Pr  -  0.69,  0.  72,  0.  75  AND  1  0 
RESPECTIVELY. 

(c)  Pr  ■  0.69,  w  =  0.67  AND  Ec  -  0.1,  0.4  AND  1  0 

RESPECTIVELY. 
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FIGURE  15.  TANGENTIAL  SLIP-VELOCITY  u^1)  IN  THE  SMALL-SLIP  LIMIT. 

(a)  Pr  =  1.0,  Ec  =  1.0  AND  co  =  0.5,  0.67,  0.  77  AND  1.0, 

RESPECTIVELY. 

(b)  Ec  =  1.0,  co  =  0.67  AND  Pr  =  0.69,  0.  72,  0.  75  AND  1.0, 

RESPECTIVE!  Y. 

(c)  Pr  =  0.69,  w  =  0.67  AND  Ec  =  0.1,  0.4  AND  1.0, 

RESPECTIVELY. 
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FIGURE  18.  NONDIMENSIONAL  SHEAR  STRESS  AT  THE  WALL  AS  A  FUNCTION  OF  THE 
NONDIMENSIONAL  DISTANCE  x. 

-  PURE-GAS  SOLUTION;  —  LARGE-SLIP  SOLUTION  FOR  THE  DUSTY 

GAS;  -  SMALL -SLIP  SOLUTION  FOR  THE  DUSTY  GAS. 
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FIGURE  19.  NONDIMENSIONAL  HEAT-TRANSFER  RATE  AT  THE  WALL  AS  A  FUNCTION  OF 
THE  NONDIMENSIONAL  DISTANCE  x. 

-  PURE-GAS  SOLUTION;  —  LARGE-SLIP  SOLUTION  FOR  THE  DUSTY 

GAS:  -  SMALL-SLIP  SOLUTION  FOR  THE  DUSTY  GAS. 


FIGURE  20.  COMPARISON  BETWEEN  THE  FLOW  PROFILES  OF  THE  GAS  WITH  AND  WITHOUT 
PARTICLES. 

(a)  VELOCITY;  (b)  TEMPERATURE. 

WITHOUT  PARTICLES;  - WITH  PARTICLES. 
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